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Line-of-sight paths are important as vhf radio channels. In a mobile 
radio system, for example, the landscape determines the communication 
possibilities in a complicated way. This paper analyzes a simple model of 
rough terrain to relate statistical terrain properties to line-of-sight paths. 
The model is constructed from conical hills, all the same heigM, distributed 
at random over the surface of a spherical earth. 

The parameters of the model are the earth's radius a, the density a 
of hills, and the grade g of the hills. Although a simpler planar model is 
obtained by letting a^<x>, a finite spherical earth is needed for most 
questions. Assuming that a base station is located at the peak of a hill, 
the most interesting line-of-sight paths are those from a typical hilltop. 
A large number of statistics of these paths are then derived, usually as 
simple functions of a, a, and g. These include properties of paths to other 
peaks, to the horizon, and to random points on the ground. 

I. INTRODUCTION 

Very-high frequency radio propagation is often said to resemble 
optical propagation. A line-of-sight path provides a good radio channel ; 
a path blocked by the terrain does not. With the aid of a topographic 
map, one can determine whether a path Q1Q2 is a line-of-sight path. 
Essentially, one must plot the ground elevation profile along the path 
to see whether the ground intersects the straight line segment QiQz. 
This calculation must include the effect of the earth's curvature. 
Atmospheric refraction is also accounted for by changing the earth's 
radius to a fictitious value. 

Having done the calculation for one path Qi, Q 2 , we learn little about 
other paths. The region covered by a transmitter at Qi, i.e., the set of 
points Q visible from Q h would be found by plotting ground elevation 
profiles along views from Qi at every possible azimuth angle. This 
region might represent the coverage of a tv station or of a base station 
in a mobile telephone system. 

This paper analyzes a statistical model to give insight into the way 
coverage regions depend on properties of the terrain. The parameters 
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of the model are the radius a of the earth, a density a of mountains 
(or hills) per unit area, and a grade (slope) g of these mountains. 
Many statistical properties of terrain and paths are then derived as 
functions of a, a, g. These properties are means, or in some cases 
distribution functions, of the random variables that appear in the 
index. Line-of-sight paths from a typical mountain peak receive special 
attention because a peak is the most likely site for a base station. 
Although the exact formulas contain integrals with unwieldy trigono- 
metric integrands, most of these formulas may be replaced by simple 
expressions, to a very good approximation. The expected area visible 
from a peak and the expected number of peaks visible from a random 
point on the ground are more complicated quantities, leading to 
integrals that are evaluated numerically. 

INDEX 

Altitude— eqs. (6), (7), (8), Table I, Fig. 6. 

Area blocking— eqs. (12), (17) to (20), (23), Table III. 

Visible— eq. (43), Table VII. 

Within horizon — eq. (35). 
Number of peaks visible : 

From a peak— eq. (26), Table V. 

From a point on the ground — eq. (37). 
Range from a peak: 

To furthest visible peak — eq. (31), Table V. 

To horizon— eqs. (33), (34), Table VI. 

To random visible peak— eqs. (25) to (29), Table IV. 

To random visible point on ground — eqs. (39) to (42), 
Fig. 15. 
Slope— Table II. 

The earth's radius a is an important parameter of the model. 
Although a simpler planar model is obtained by letting a — > °o f the 
planar model is inadequate for most statistics of interest. 

With a and a fixed, the terrain becomes rougher as g increases. As a 
rule, the model predicts more long line-of-sight paths and larger 
expected visible area for rougher terrains. However, in mobile radio 
these long paths are more important as sources of interference than 
as useful channels. 

II. THE MODEL 

The terrain model will use conical mountains distributed at random 
in a Poisson pattern over the surface of the earth. Begin with a sphere 
of radius a miles (o may be the true radius of the earth, or something 
larger if atmospheric refraction effects are to be taken into account). 
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Place points at random on the surface S of this sphere using a Poisson 
process with density a points per square mile. Each Poisson point will 
represent a mountain peak, and so the sphere of radius a will be called 
the peak sphere. 

Each Poisson point P will be associated with a mountain-shaped 
subset M(P) of the interior of the peak sphere. The subsets M(P{), 
M(P 2 ), ■ • • for the various peaks will overlap. Take the union of all 
the subsets M{P) to represent the earth. 

The simplest shape for M{P) is the cone consisting of all rays from 
P making angle <0 with the inward-pointing normal to S. This cone 
has to be truncated to keep it from extending beyond the peak sphere 
in the direction antipodal to P. The surface of the cone is tangent to 
an inner sphere, concentric with the peak sphere and having radius 
a sin 0. Take M (P) to be the inner sphere plus the part of the cone 
that lies between P and the inner sphere. Figure 1 shows M(P) 
shaded. 

With this construction, the terrain consists of conical mountains, 
all having the same height and the same grade g = cot 0. There may 
also be flat places where the earth's surface coincides with the inner 
sphere. A flat spot occurs at any point that lies further than (§*• — 0) 
radians away from all Poisson points. Flat spots are rare, except when 
the parameters a, are chosen to produce widely separated mountains 
having very gentle slope. 




ANTIPODE 



Fig. 1 — Construction of M(P). 
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Figure 2 is an elevation contour map for a typical random terrain. 
Some unrealistic features of the model are evident. The conically 
shaped mountains have circular contour lines. The peaks are distrib- 
uted chaotically instead of being arranged in rows (mountain ranges). 
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Fig. 2 — Contour map. The symbols -8+6 — 4:2 denote altitude levels ordered 
from the peak sphere downward. 
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Figure 3 is a plane cross section through the earth in the same model. 
This figure is more interesting for the present problem because the 
existence of a clear line-of-sight path between two points depends only 
on the shape of such a profile. Note that the elevation curve in Fig. 3 is 
composed of convex arcs (hyperbolas) that join in the valleys between 
mountains. But, at least, the maxima in Fig. 3 have different heights. 
Figure 4 shows random terrain as seen from one of the peaks looking 
out toward the horizon. The nearest and furthest peaks shown have 
ranges of 6 and 150 miles. The parameters were picked to match a 
particular portion of the Alps for which a panoramic photograph was 
also available. The deficiencies of the model are less evident in this 
figure. The curvature of the earth makes it less obvious that all peaks 
have the same height. 

In real terrain, it is sometimes possible to see part of a mountain 
even though the mountain's peak is obscured from view. That cannot 
happen in this model, as will now be proved. Suppose that the view 
of a peak Pi is blocked when the eye is at E. Then the line segment 
P\E contains a blocking point B 2 belonging to another mountain 
M(P 2 ). Now consider any other point P of M(Pi). P must lie on some 
line segment PJ, where / belongs to the inner sphere. Figure 5 shows 
the triangle EPJ. The segments EP and B 2 I cross at some point B 
in the triangle. B belongs to the convex set M(P 2 ) because B 2 and / 
belong. Then B is a point of M(P 2 ) blocking the view of P. 

By making a — * °° , one obtains a planar model of random topography. 
The peak sphere S becomes a peak plane. At a point Q, the land surface 
lies below the peak plane a distance 

y = g Min ||P« - Q\\, (1) 

i 

where the minimization is over all Poisson points P,. Replacing S by a 
plane simplifies the analysis considerably but, unfortunately, it 
produces a much less realistic model. If Fig. 4 has been drawn for a 
planar model, every peak P, would have been visible. Even worse, 
Section VIII shows that the expected area visible from a peak would 
be infinite. For that reason, the extra complication of a spherical earth 
is really necessary for some questions about line-of-sight paths. 

III. PARAMETER ESTIMATION 

The two parameters <r, g = cot 6 can be chosen to fit the model to 
terrain measurements. One might estimate the density a by counting 
peaks. A difficulty is that one must then decide how big a hill must be 
to be counted. Surely every bump on the landscape ought not to count 
as a peak. This decision is avoided by using statistical properties of 
the point Q lying below a random point q on the peak sphere. The 
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I 
Fig. 5— M(P 2 ) blocks all of M (Pi) from view if it blocks Pi. 

altitude and slope of the terrain at Q are two useful random variables. 
Both depend on the angle 7 = < POQ from Q to the nearest peak P 
(see Fig. 1). Since a circular cap of angle 70 on the peak sphere has 
area 2ira 2 (l - cos 70), the distribution function for <POQ can be 
written immediately, 

Prob (7 ^ 70) = 1 ~ exp [- 2wa t a(l - cos 70)]- (2) 

There is no natural sea level in the model, and so it will be con- 
venient to specify the altitude at Q by giving the depth y, measured 
from 8 down to the land. If 7 ^ fa - 0, then ground level coincides 
with the inner sphere, i.e., 

y = a(l - sin 0), y ^ fa - 6. (3) 

For smaller angles 7, 

y - all- sin 0/sin (7 + 6)2, £ y < fa - 6, (4) 

as is clear from Fig. 1. These formulas, together with the distribution 
(2) for 7, determine the depth distribution, 

Prob [y £ a[l - sin 0/sin (7 + 0)]} 

= 1 - exp [-2*aV(l - cos 7)], Ogy <fa-8 (5) 
Prob [y ^ a(l - sin0)} = 1. 

Although one can easily tabulate the distribution function for y by 
substituting numerical values of 7 into (5), the distribution function 
is easier to visualize in a limiting case. Since a is a large radius, let 
a— +00 in (5). As one might expect, the formulas tend toward the 
depth distribution function in the planar model, 

Prob {y ^ Y] = 1 - exp \_-<nr{Y/gY2- (6) 

In this limit, a and g enter the distribution only via a single length 
parameter a^g, which is an index of altitude variability. Thus, altitude 
distribution data alone cannot be expected to supply good estimates of 
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both g and a. Some simpler statistics are the median, 

Median (y) = (tt" 1 log e 2g 2 /a)i 
= 0.4697<r-ty, 

and the moments, 

E(y") = r(l + *fc)(<77™)*' 2 . 
Particularly, the mean is 

y = E(y) = \^g (7) 

and the standard deviation is 

[Var (y)y = [(tt" 1 - 2-»)0«/cr]» = 0.2683«7-ty. (8) 

It is also possible to obtain (6) as an exact result for a spherical 
model in which the shape of the mountains is only approximately 
conical. That entails a new choice of the set M(P) in Fig. 1. Define the 
new shape so that the depth becomes 



y = 2ga sin §7, 



(9) 



where again 7 is the angle to the peak. At P, M(P) comes to a point 
approximating a cone of slope g. At the antipode to P, M(P) has depth 
2ga ; then this model requires g < \. Now the depth distribution for all 
7 is again given by (5) but with the left-hand side replaced by 

Prob \y ^ 2gasin^y}. But that is (6), exactly. 

For many values of g and a, the planar approximation (6) to the 
depth distribution (5) is very good. For example, Table I compares 
the planar approximation with some distributions having a = 3959 mi, 
the earth's radius. In the table, the cones have grades g = 0.05, 0.1, 
and 0.2 and the density a is adjusted to fix the standard deviation in 
(8) at 528 ft (0.1 mi). Table I gives percentiles of the distribution as 

Table I — Altitude percentiles (in feet) 



Spherical Model 


Planar 


9 = 
a = 


0.05 
0.0171 


0.1 
0.0683 


0.2 
0.2732 


o/g i = 6.831 


0.1% 


-1899 


-1964 


-1980 


-1986 


1% 


-1378 


-1421 


-1432 


-1436 


10% 


-691 


-712 


-717 


-719 


25% 


-315 


-327 


-331 


-332 


50% 


70 


63 


62 


61 


75% 


402 


400 


399 


399 


90% 


642 


641 


640 


640 


99% 


896 


896 


896 


896 
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altitudes measured upward from a common level, corresponding to the 
depth y in (7). 

Figure 6 is an altitude distribution for northern New Jersey. It 
was obtained from a topographic map by reading altitudes at 52 
points, 10 km apart in a rectangular grid covering latitudes 40°30' 
to 41° and longitudes west of 74°. The altitudes ranged from to 
1100 ft. Data for parts of New Jersey further south were not used; the 
topography of New Jersey is too variable for both north and south 
to be well represented by a single simple model. The planar model 
fits the observed points well, except at low altitudes. As an alternative, 
use the spherical model with a = 3959 mi. By taking g = 0.011, one 
obtains a maximum depth (3) near 1100 ft, so that low altitudes can 
be regarded as occurring on the inner sphere. Then a remains as a 
parameter to adjust for a good fit. 

The parameter g = cot is the grade at mountain peaks. At the 
random point Q, at angle y away from a peak in Fig. 1, the grade is 
smaller because the normal to the conical surface makes an angle 
\ir — d — y with the vertical direction OQ. Thus, the grade at Q is 

, _ jcot (0 + 7) - (9 - tan 7 )/(l + g tan 7 ), if B + y ^ \w 
[0 otherwise. 




100 200 300 400 500 600 700 800 900 1000 1100 
ALTITUDE IN FEET ABOVE SEA LEVEL 

Fig. 6 — Altitude distribution for northern New Jersey. Curve is for planar model 
with peak sphere at 1130-ft altitude and [Var (j/)]» = 400 ft. 
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Table II — Percentiles of G/g 



small 0.1 0.2 0.5 



0% 


-1.0000 


-1.0000 


-1.0000 


-1.0000 


1% 


-0.9995 


-0.9995 


-0.9995 


-0.9994 


10% 


-0.9510 


-0.9506 


-0.9492 


-0.9399 


25% 


-0.7070 


-0.7053 


-0.7001 


-0.6667 


50% 














75% 


0.7070 


0.7053 


0.7001 


0.6667 


90% 


0.9510 


0.9506 


0.9492 


0.9399 


99% 


0.9995 


0.9995 


0.9995 


0.9994 


100% 


1.0000 


1.0000 


1.0000 


1.0000 



The grade occurs on the inner sphere. This result, together with (2), 
determines the distribution of the grade g'. In most cases, the grade g' 
has high probability of being close to g; one should not expect this 
distribution to fit observed grade data well. 

At q, one might move in a random direction and ask for the slope 
along the random path through Q. The slope, which depends on the 
angle <p between the path direction and the uphill direction, lies in the 
range — g'^G^g'. With some simple geometry, one finds 

G = g' cos <p/(l + /2 sin 2 *>)*. 

By using the known distribution for g' and assuming a flat distribution 
for <p, one can obtain a distribution function for G. This would be the 
distribution of the slopes G seen in cross sections like Fig. 3. A simple 
distribution is obtained only in the planar model limit, for which 
g' = g = cot identically : 

Prob [G ^ tan X} = 1 — ir -1 arc cos {sin x/cos 6\. 

Table II gives the slope distribution in the planar model for several 
values of g. In the limit of small g, the distribution function for G/g 
tends to 1 — tt -1 arc cos G/g. 

IV. BLOCKING REGIONS 

Suppose two points Q h Q 2 are given, representing the positions of 
two antennas. In general, Q h Q 2 can He anywhere above the inner 
sphere. A clear line-of-sight path exists between Q x and Q 2 as long as 
the straight-line segment Q1Q2 does not intersect any of the sets 
M(Pi). The blocking region for Q h Q 2 is the (open) set of points P on S 
such that QiQ 2 intersects M(P). The area of the blocking region enters 
into the probability that a line-of-sight path Q\Q 2 exists. The advantage 
of the conical mountains M(P) is that blocking regions assume simple 
shapes. 
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The simplest blocking region is one for a pair of points Q lt Q 2 both 
on S. If the line Q1Q2 intersects the inner sphere, all M{P { ) block the 
path. The blocking region then consists of the entire sphere S. If 
Q1Q2 misses the inner sphere, then blocking occurs at a point Q on the 
path Q1Q2 if a peak P, lies too close to Q. If the depth of Q is y, then 
(4) gives the angle 7 to peaks P such that Q lies on the surface of 
M(P). Then blocking occurs at Q if a circular cap of angular radius 
7 contains a peak. The pole of this cap is the radial projection q of Q 
onto S. The blocking region for the path Q1Q2 is the union of all the 
blocking caps for points Q on the path. These caps are largest midway 
between Qi and Q 2 , shrinking to points at Qi and Q 2 . Then the blocking 
region is lens-shaped, as in Fig. 7. 

Figure 7 shows two arcs K, K' which form the boundary of the 
blocking region. The argument that follows shows that K, K' are 
actually arcs of circles. Figure 8 is another view of the peak sphere 
projected directly along the line Qi, Qi. Two planes, ir and ir', can be 
drawn through Q h Q 2 and tangent to the inner sphere, say at C and C". 
These planes project to lines in Fig. 8. The planes ir and ir' intersect S 
in two circles, centered at C and C and both passing through Qi and 
Qz. Since M(P) is the convex hull of P and the inner sphere, ir is a 
supporting plane of M{P) as long as P lies below t (i.e., in the half- 
space containing the inner sphere). Then M(P) does not block the 
path QxQi if P lies below ir, or below v'. The part of S lying above 
both ir and ir' appears shaded in Fig. 8. Suppose P belongs to the shaded 
region. Project the triangle C'PC, a subset of M(P), onto the plane 
of Fig. 8. The path QiQi projects to a point lying inside this projected 
triangle. Then Q1Q2, a chord of S, must intersect the triangle C'PC. 




Fig. 7 — Blocking region for two points Qi, Qt on the peak sphere. 
746 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1975 




Fig. 8 — Another view of the blocking region. 

The point of intersection is a point of M(P), which blocks the path. 
Thus, the shaded region, bounded by arcs of the circles S (~) t, S f) ir', 
is the blocking region for Q1Q2. 

The area A (Qi, Q 2 ) of the blocking region in Fig. 7 will now be 
expressed as a function of the angle 2p = zQ\OQi. Project the centers 
C, C in Fig. 8 radially out to c, c' on S. Figure 9 is another view of S 
showing c, c' as the poles of two circular caps bounded by S f) «" an d 
S (~) t'. The angular radius of both caps is fir — 6, as is clear from 
Fig. 8. The chord Q1Q2 subtends some angle 2a = z QicQz at c. Using 
the spherical sine law in the right triangle Q\, c, |(Qi 4- Q2), one may 
determine a from 

sin a = sin p /cos 6. (10) 

The cap with pole c has area 27ra 2 (l — sin 6) and the sector included 
within angle 2a has area 

A a = 2aa 2 (l - sin0). 

Also, the triangle QicQz has area 

A T = (2a + 2/3 - ir)a 2 , 

where /3 = Q1Q2C1 = QzQic. The sine law may be applied to triangle 
QicQ* to find /3 

sin /3 = cos 6 sin 2a/sin 2p. (11) 
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Fig. 9 — Angles used in deriving area A of blocking region. 

The difference A s - A T is hMQh Qa). Thus, 

A(Q h Q*) = a 2 (27r - 4/3 - 4a sin 6), 



(12) 



where (10), (11) give a, 0. 

The blocking region is more complicated if Q\, Q 2 or both are not 
on S. As in Fig. 7, each point Q on Q1Q2 is blocked by peaks lying in a 
circular cap of radius 7 given by (4) ; the blocking region is the union 
of these caps. Let Qi, Q 2 be the points where the extended line Q1Q2 
meets S. The blocking region for Q1Q2 is a subset of the blocking 
region for Q[Qi As shown in Fig. 10, the blocking region consists of 
the caps for blocking at Qr and Q 2 plus the part of the blocking region 
for Q1Q2 that lies between these caps. The centers of the two end caps 
are the points q h q 2 obtained by projecting Q h Q 2 radially onto £. 

The two end caps have a special role in the blocking. Normally, 
Qi, Q2 are known to lie above ground, and so the two end caps are 
known to contain no peaks. If the ground levels below qi, q 2 are known, 
then peaks P lf P 2 must exist somewhere at the appropriate angles 71, 
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72 away from q h q 2 . In Fig. 10, Q h Q 2 are assumed to be at ground 
level ; then Pi, P 2 lie on the boundaries of the end caps. The mountains 
M(Pi), M(P t ) on which Q h Q 2 lie can themselves block the path QiQ 2 . 
Thus, in Fig. 10, P 2 blocks the path because it lies in the blocking 
region. To compute the conditional blocking probability for the 
configuration in Fig. 10, one must know both the area of the part of 
the blocking region that lies outside the end caps and also angles <pi, <p 2 
that limit where Pi, P 2 can lie to cause blocking. In the applications 
that follow, it will suffice to let Qx lie at a peak Qi = Pi and let Q 2 be a 
point at ground level. That simplifies Fig. 10 to Fig. 11. 

Let z 2 = z Pi0q 2 , the angular distance along the arc Piq 2 . The depth 
at Q 2 determines the angle 72 of the end cap. The sides of the spherical 
triangle P\q 2 c are now known, and so its angles /3 = LqiPic, 
f2 = ZqtfPi, ir — <p 2 = /_P\q 2 c are determined. One finds 



tan \ tp 2 = 



1 — cos (z 2 — 72) + g sin (g 2 — 72) 
g sin (22 + 72) — 1 + cos (z 2 + 72) 

sin jQ = sin <p 2 cos (6 + 72) /cos 6 

sin f 2 = sin ip 2 sin z 2 /cos 6. 



(13) 

(14) 
(15) 



The blocking area is twice the area of the half of the blocking region 
above the line P\Q' 2 in Fig. 11. That half can be obtained as a sum of 
two parts. One part is a sector of angle tr — <p 2 from the end cap ; its 
area is (w — <p 2 ) (1 — cos 72) a 2 . The other part is obtained by removing 




Fig. 10 — Blocking region for points Qi, Qt which are not peaks. 
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Fig. 11 — Blocking region for Qi = P, a peak, Q2 not a peak. 

the triangle Piqic of area (/3 + f 2 — ^2) a 2 from a sector centered at c. 
The sector has area f 2 (l — sin 5) a 2 . These areas may be combined to 
express the area of the blocking region in the form 



where 
and 



A(Qi,Q 2 ) = A Q + A 2 , 
Ai = 2tt(1 — cos 7 2 )a 2 



(16) 



A = (2^2 cos 72 - 2/3 - 2r 2 sin 9)a 2 . (17) 

A 2 is the area of the end cap and A is the area of the remainder of the 
blocking region. 
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Although the blocking area A (Q h Q 2 ) for the general situation of 
Fig. 9 will not be needed, it can be obtained in the form 

A(Qi, Q a ) = A(Q[, Q t ) + A(Q h Q 2 ) - A(Q' h Q' 2 ). (18) 

Note that formulas like (16) and (17) give A(Q' h Q 2 ), A(Q h Q 2 ) while 
(10), (11), and (12) give A(Q[, Q' 2 ). Likewise, with a change of sub- 
scripts, (13) gives <pi as well as <£> 2 . 

These formulas can now be used to obtain the path probability 
p(Qh Qz), the conditional probability that a clear line-of-sight path 
exists between given points Q h Q 2 . When Qi, Q 2 are on S, as in Fig. 7, 
V (Qh Q2) is J ust the probability that the shaded region of area A (Q h Q 2 ) 
contains no peaks. Then 

p(Qi, Q2) = exp l-<rA(Qi, Q,)], (19) 

with A(Qi, Qt) given by (12). 

The situation in Fig. 11 is more complicated. Qi = Pi, a peak, and 
Q 2 is supposed to lie on the ground. Then the cap of area A 2 is known 
to be empty. Two conditions must hold if the entire blocking region 
is to be empty. One is that the peak P 2 of the mountain on which Q 2 
lies causes no blocking. Since P 2 is equally likely to be anywhere on 
the boundary of the cap around Q 2 , there is probability 1 — <p 2 /ir that 
P 2 does not block the path. The second condition is that the remainder 
of area A of the blocking region is empty. Then 

p(Qh Qi) = (1 - vtlr) exp (-aA ), (20) 

where (13) and (17) give <p 2 and A . Formula (20) applies as long as 
72 < z 2 . It is also possible to have y 2 = z 2 . In that case, Q 2 lies on the 
mountain M(Pi) ; the path in question runs from the peak Q\ = P\ 
to Q 2 along the surface of the cone M(Pi). Whether or not such a path 
is to be considered blocked is a matter of definition. Here Af (Pi) is 
regarded as an open set so that the path is not blocked. As 72 — * z 2 , 
one finds <p 2 — >0 and A — >0 so that p(Q\,Qi) — > 1, i.e., (20) con- 
tinues to give the correct probability in the limit. 

Another limiting situation, 72 — > 0, illustrates an important distinc- 
tion between Figs. 11 and 7. In the limit A 2 — > and A becomes the 
area of a lens-shaped region, such as shown in Fig. 7. Then the ex- 
ponential factor in (20) becomes the path probability (19) for the 
two peaks Qi( = Pi), and lim Q 2 . However, (20) contains an extra 
factor (1 — (p 2 /ir) which approaches \, not 1. This disagreement 
between (19) and (20) is explained as follows. From a point Q 2 near a 
peak P 2 , one can look over a 180-degree view; the mountain M(P 2 ) 
blocks the other 180 degrees. Then the factor \ in (20) is needed to 
account for possible blocking by M{P 2 ). But exactly at the peak 
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Q„ = p 2) the mountain M(P 2 ) no longer interferes in any direction. 
Then no factor \ is needed in (19). The discontinuity in p(Q h Q 2 ) as 
q 2 _> p 2 could be avoided by assuming that the antenna location Q 2 
lies at some known positive height h above ground. 



V. PATHS BETWEEN PEAKS 

The simplest blocking regions were for paths Q1Q2 with both end 
points on mountain peaks. The path probability p(Q h Q 2 ) in (19) can 
now be used to derive some interesting properties of peak-to-peak 
paths. In this section, Qi will be a given peak Pi. Q 2 will be another 
peak selected at random. An element of area dA(Q 2 ) on the peak 
sphere S has probability adA(Q 2 ) of containing a peak Q 2 . Then 
0p(Qii Qi)dA(Qi) is the probability that the element contains a peak 
Q 2 which is visible from Q\. 

Let d{Qi, Q 2 ) denote great circle distance between Qi and Q 2 . Let 
2*(d) denote the random variable which is the sum 



Md) - L d(Q h Pi) 1 



(21) 



of fcth powers of distances from Qi to all other visible peaks P< lying 
within distance d[d{Q h Pi) ^ d~\. The element dA(Q 2 ) contributes a 
term d(Q h Qt) k to S*(d) with probability trp(Qi, Q*)dA(Q 2 ). Thus, the 
expected value of 2 fc (d) is 

ElMd)l = * J J d(Q h Q z ) k p(Qi, Q2)dA(Q 2 ), (22) 

where the integral extends over all points Q 2 in the cap d(Q\, Q 2 ) ^ d. 
Another random variable Z k is a sum like (21) extended over all 
visible peaks, at any distance from Q x . The mean E(2 k ) is an integral 
(22) over the entire sphere. Evaluating (22) will give the mean number 



Table III — Blocking area 4(Qi, Q2) in square miles, as given by 

exact and approximate formulas (12) and (23) with 

range d{Q lt Q 2 ) = 100 miles 



Grade g 


Exact 


Approximate 


0.01263 


7887.7 


3333.1 


0.02 


2430.9 


2104.9 


0.05 


857.6 


842.0 


0.1 


423.2 


421.0 


0.2 


210.2 


210.5 


0.5 


84.1 


84.2 


1.0 


42.0 


42.1 
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of visible peaks 2? (So) and other information about the distances to 
visible peaks. That could be done numerically, using (10), (11), (12), 
and (19). However, the approximation that follows simplifies the 
evaluation. 

The approximation is one which holds when a is so large that angles 
2p between visible peaks can be considered small. The planar model 
has A(Q h Qi) = and p(Qi, Qi) = 1, which is too rough to make 
sense in (22). Instead, the first nonzero term in a series for A(Qi, Qi) 
in powers of p will be used. Expansion of the exact formulas (10), (11), 
and (12) is laborious but straightforward: 

« = p(i + g 2 ) h /g + p 3 (i + g 2 W(W) + O(p') 

= J* - P /g - W + 1)p7(6<7 3 ) + 0(p 8 ) 
A(Q 1( Q 2 ) = a 2 (2p)7(60) + 0(p 5 ). 
= r*/(Gga) + ■■•, 

where r = 2ap is the great circle distance d(Q h Q 2 ). 

For a simpler, more intuitive, derivation, one may find the size of 
the circle about a typical point q along the path Q h Q 2 in Fig. 7. If z is 
the great circle distance from Q\io q, then the ground level below q 
lies at depth y satisfying 

(a — y) cos (z — \r)/a = a cos \r/a 

or 

y = z(r- g) /(2a) + •••. 

The radius ay of the circle about q is approximately y/g, and the 
blocking area is approximately 

A(Qi, Qi) = J\ydz 

= [' z(r- z)dz/(2ag) 
Jo 

A(Qi,Q 2 ) =r«/(60a)+ ■■-, (23) 

as before. 

From the form of the series used in deriving (23), one may predict 
that the rate of convergence is determined by the ratio p/g. Table III 
shows that (23) does give a better approximation for large g than for 
small. In Table III, a = 3959 mi; a large range, 100 mi, was used for 
a severe test of the approximation. At grades g smaller than 0.01263, 
a 100-mi path between peaks is blocked by the inner sphere. The small 
p/g condition is another way of requiring that the path QiQ 2 clears 
the inner sphere by a safe margin. 
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Now use the approximation (23) for A(Q h Q 2 ) in (19) to evaluate 
the integral (22) for the expectation [2*(d) ] : 

rdlita) 

tf[2*(d)] = 2™ / (2a P )* +1 exp { - 4a 2 ap 3 /3? } 2adp 

Jo 

ELMd)~] = (2tt(7/3)D* +2 / «(*-»)/• exp (-«)<**, (24) 

where 

Z) = (6a<7/<r)*. 

The integral in (24) may be expressed in terms of the incomplete 
gamma function, 

E\J, k (d)~\ = (2W3)£ fc+2 {r[(/c + 2)/3] - r[(fc + 2)/3], (d/£>) 3 }, 

or the X 2 distribution function, 

E&t(d)l = (2W3)i) fc+2 r[(fc + 2)/3]P(X»|v), 

where 

x 2 = 2(d/D) 3 

and the number of degrees of freedom is 

v = 2(fc + 2)/3. 

Although the approximation (23) becomes poor at long ranges, the 
integrand is very small there. Thus, (24) can be expected to hold even 
for long ranges. In particular, the expectation E(2 k ), for visible peaks 
of all ranges, may be approximated by letting d — ►«> in (24) : 

E(2 k ) = (2x«7/3)D t+2 r[(fc + 2)/3]. (25) 

For the special value k = 0, (25) gives the mean number of peaks 
visible from Qi : 

E(2 ) = 7r(rZ) 2 r(5/3) 

= 9.3645 (aV<r)* 

= 2344 (g 2 <r)* if a = 3959 mi. (26) 

Note, as predicted earlier, that the mean number of visible peaks 
tends to infinity in the limit of large a (planar model). When k = 1, 
(25) simplifies to 

E(2i) = Irag. (27) 

As a increases, (26) shows that the mean number of visible peaks 
increases, but (27) shows that the mean sum of distances to visible 
peaks remains unchanged. This indicates that visible peaks tend to 
be closer for large a than for small. One way to define a range for a 
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"typical" peak is to form the ratio 

Dx = B(2d/B(2J = D/r(2/3) 
= 0.738487D 
- 1.3419O(o0/«r)*. (28) 

VI. RANGES BETWEEN PEAKS 

One might ask for a probability distribution for the range d from 
a peak Pi to a randomly chosen visible peak P ^ P\. The random 
process for choosing a peak must be specified with care. Perhaps the 
most natural process would be this. Construct a random landscape 
and choose a peak P from the set of 2 visible peaks, all peaks equally 
likely. Then ask for the probability that P is one of the 2 (d) peaks 
within range d of Pi. Given a landscape, the conditional probability 
that P is within range is 2 (d)/2 . Then the unconditioned probability 
is £ , [2(d)/2 ]. Unfortunately, the expectation is hard to obtain [there 
is also a question of giving an appropriate meaning to 2 (rf)/2 when 
2„(d) = 2o = 0]. 

By using a different random process, one obtains a simpler distri- 
bution. Construct a trial random landscape and pick one of the peaks 
P at random, this time from the set of all peaks on the entire sphere S. 
P may not be visible. If not, discard that trial and construct a new 
landscape. Continue constructing landscapes and choosing peaks until 
the chosen point P is visible. Then ask for the probability p(d) that 
P lies within range d. 

To determine p(d), note that the total number of peaks on the entire 
sphere has the Poisson distribution with mean 4fira 2 <r. The argument to 
follow assumes that this number is large, so that the number of peaks 
actually obtained is almost always very close to its mean value. Then 
the probability that P is visible is £(2 )/(47raV). If g[2 , 2 (d)] is the 
joint probability for 2 and 2 (<2) in each trial, then <?[2 , 2 (d)]2 / 
(4ira 2 a) is the probability that a trial has 2 visible peaks, 2 (d) within 
range d, and that P is visible. The joint probability for the numbers 2 , 
2 (rf) of the landscape, selected when P is visible, is g'[ 2 o, 2 (d)] 
= g[2 , 2 (rf)]2 /^(2o). The probability that P lies within range d 
is obtained as a sum over 2 (d) and 2 

P(d) = L g'[2o, 2 (d)I|2o(d)/2o 

= £[2o(d)]/£(2 ) 
p(d) = 1 - T[2/3, (d/.D) 3 ]/r(2/3), (29) 

the last line following from (24). 

It is clear from this derivation that the second random process 
tends to select random landscapes with larger 2 than the landscapes 
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Table IV — Probability p[d) = E(2 (d))/E(2o) that a randomly 
chosen visible peak lies within range d 



d/D 


Probability 


0.25 


0.06880 


0.30211 


0.1 


0.48595 


0.25 


0.5 


0.26361 


0.72212 


0.5 


0.75 


0.53050 


1 


0.77518 


1.20507 


0.9 


1.32182 


0.95 


1.5 


0.98440 


1.55886 


0.99 


1.81350 


0.999 


2 


0.99983 



of the first process. However, S may be expected to have a highly 
peaked distribution, in which case 2 is nearly always close to E(Z ). 
Then q(-, •) and q'{-, •) are nearly the same, and (29) is also a good 
approximation to the range distribution for the first random process. 

Equation (29) provides numerical values for the range distribution 
in Table IV. 

Another random variable of interest is the range to the furthest 
visible peak. Even the expectation of this maximum range is hard to 
derive. However, a simpler "typical" maximum range is the range d m 
such that the expected number of visible peaks with ranges d > d m is 
just \. Then d m satisfies 



tf(Zo) - E[2o(d m )l = h 



(30) 



and (24) shows that 



/ 



(dm/D) 1 



*exp(-w)du = 3/(4tt(tD 2 ). 



(31) 



Table V — Mean number of visible peaks E(2 ) and range d, 
such that E(2 -2 (d m )) = V2 



<tD* 


a 2 g 2 ir 


£(2o) 


d m /D 


3.11 


0.84 


8.8 


1.30 


5.69 


5.12 


16.1 


1.40 


11.3 


40.1 


32.1 


1.50 


24.5 


409 


69.5 


1.60 


58.4 


5528 


165.6 


1.70 


153.6 


100572 


435.5 


1.80 


450.5 


2.54 X 10 6 


1278 


1.90 


1477 


8.95 X 10 7 


4189 


2.00 


5447 


4.49 X 10» 


15448 


2.10 
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Table V gives numerical values of d m /D as a function of <tD 2 
= (36o 8 flfV)*- #( s o), which also depends on <tD 2 as shown by (26), 
also appears in the table. Note that d m is not just a function of a single 
product of powers of a, g, a; it has a more complicated form (ag/a)* 
X function (a 2 g 2 <r). 

The integral in (31) is a rapidly decreasing function of d m /D. Then 
the numbers in Table V would not change much if d m were redefined 
with the term § in (30) replaced by any other number of the same order 
of magnitude. For the same reason, d m can be expected to be a good 
approximation to the mean range to the furthest peak. 

VII. THE HORIZON 

The approximation (23) will now be used to derive properties of 
the range from a peak Pi to the horizon at a random azimuth angle. 
The range to the horizon is a more interesting random variable than 
the range to a random visible peak. As has been noted, it is not always 
clear what to count as a peak in a real landscape. But the horizon 
has no ambiguity. 

Look from Pi with a fixed azimuth angle. One sees only sky at high 
elevations and ground at low elevations. The horizon point is the limit- 
ing point at ground level which has the highest elevation angle. The 
distance z from Pi to the horizon is the range of interest here. 

Figure 11 will now be used to derive the conditions under which Q 2 
is the horizon point, as seen from a peak Pi. If Q 2 is the horizon point, 
the entire straight line path PiQ 2 in Fig. 11 must intersect the ground 
only at Q 2 . Then the entire lens-shaped blocking region for Q 2 must 
contain no peaks. But the depth at Q 2 determines the circle on which 
a peak must lie. This circle appears in Fig. 11 inside the (open) blocking 
region. The only place that this peak can be now is on the boundary 
of the blocking region at one of the points of tangency T, T'. 

Figure 11 shows the usual situation in which the horizon point is 
not on the inner sphere. There is small probability that Q 2 is on the 
inner sphere. In that case, q 2 is at angle %tt — 6 away from Pi, the 
centers c, c' coincide with q 2 , and the blocking region is bounded by the 
circle through Pi with center q 2 . There is no second peak on the bound- 
ary of the blocking region; Qi lies on M (Pi). 

To find the probability distribution for the horizon range, one may 
first find the joint distribution for that range and the range r to the 
intersection point Q 2 . Since r is the largest range for which the corre- 
sponding blocking region is empty, the probability distribution func- 
tion for r is P(r) = 1 — exp [— a A (Pi, Q' 2 )~\. In this derivation, the 
approximation (23) will be used to write 

P(r) = 1 - exp [- (r/D) 3 ]. 
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Figure 12 shows Q' 2 lying at a range between r and r + dr, an event 
of probability dP(r). Given this position for Q 2 ', the band between the 
boundaries of the blocking regions at r and r + dr contains the peak 
on which the horizon point lies. The shaded part of this band is the 
region in which the peak must lie so that the horizon point Q2 will 
have range 2 2 ^ z. The conditional probability function for the horizon 
range is just the ratio of the area of the shaded part of the band to the 
total area of the band. To simplify that calculation, one may replace 
the dotted line by a great circle that crosses P\P' 2 at right angles. That 
approximation leads ultimately to the conditional distribution 



Prob {horizon range ^ z\r) = (z/r) 2 , ^ z ^ r. 



(32) 



The details are omitted because the result can almost be guessed 
immediately from the roughly triangular shapes of the two parts of the 
shaded region. 

Now the unconditional probability distribution for the horizon range 
is obtainable from (32) by integrating 

Prob {horizon range ^ 2} 

= f* dP(r) + J* (2/r)*dP(r) 

= (z/D)*m, (z/Dyi + 1 - exp [-(2/D) 3 ]. (33) 
Table VI gives numerical values. 




r + d r 



Fig. 12 — Horizon at range ^ 2. 
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Table VI — Distribution function for range z to the horizon 
looking from peak P x with a random azimuth angle 



z/D 


Probability 


0.21417 


0.1 


0.25 


0.13625 


0.35618 


0.25 


0.5 


0.42355 


0.56305 


0.5 


0.75 


0.70432 


0.79977 


0.75 


1.0 


0.88853 


1.02324 


0.9 


1.15749 


0.95 


1.25 


0.97109 


1.40527 


0.99 


1.5 


0.995247 


1.67110 


0.999 



The moments of the horizon range z are easy to find. From (33), the 
probability density for z is 

2z r r~HP{r). 

The fcth moment of z is 

E{z k ) = 2 P z k+1 f" r- 2 dP(r)dz 
= [2/(* + 2)]£(r*). 

The last line is obtained by integrating by parts. The expectation on 
the right is another integral that can be evaluated in the manner of 
(24) and (25). The final result is 

E(z k ) = 2D*r[l + (fc/3)]/(fc + 2). (34) 

Equation (34) with A; = 2 is particularly interesting. If z(<p) is the 
range to the horizon when looking with azimuth <p, then the mean 
area within horizon range is 

E (area within horizon) = Ey- I z 2 (<p)d<p) 

= icE{z 2 ) 

= §xr(5/3)D 2 

= 1.41803Z) 2 . (35) 

This expectation is only an upper bound on the mean area visible. 
For, as is clear in Fig. 4, there are points within horizon range that are 
obscured from view. 
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It is interesting to compare Tables IV and VI. At any given proba- 
bility level, the range to the horizon is smaller than the range to a 
randomly chosen visible peak. This may be surprising at first. How- 
ever, each visible peak is itself a point on the horizon. As seen in Fig. 4, 
the horizon consists entirely of small line segments extending down the 
sides of the mountains from the visible peaks. The line segments for 
distant peaks tend to subtend smaller azimuth angles than the seg- 
ments for nearby peaks. Picking an azimuth at random, one is more 
likely to find the horizon point on one of the nearby visible peaks 
than on one far away. 

Another expectation that exhibits the same effect is the mean 
azimuth angle between the horizon point and the peak of the mountain 
on which the horizon point lies. The ranges z and r determine this 
angle. Without belaboring the details, one can approximate this angle 
by its tangent and make the further approximations by which (33) 
was derived. The expected angle is found to be 

£(angle) = E[(r - z)/(2ag)] 
= T(4/3)D/(6ag). 

That result can be stated in a more illuminating way : 

^(ZoWangle) = *rr(4/3)r(5/3) 
= (2tt3-*) (2?r) 
= 0.40306 (2*-). 

By contrast, if E(2 ) peaks were evenly distributed in azimuth with 
angular separation 2ir/E(Z ), one would obtain 

E (So) E (angle) = 0.25 (2tt). 

The larger factor 0.40306 again occurs because of the variability of the 
angles which visible mountains subtend on the horizon. 

VIII. COVERAGE AREA 

The coverage set for a point P is the set of points Q such that a line-of- 
sight path PQ exists. In vhf radio applications, the coverage set of P 
is the set of points Q to which an antenna at P can radiate a strong 
signal. This section will estimate the mean coverage area C, the expected 
area of the coverage set for P = Pi, a peak. 

The coverage set can have a very complicated shape. Figure 13 
shows one coverage set. In Fig. 13, the peaks are not in a Poisson 
pattern ; to simplify the drawing, the peaks were located at points of a 
regular lattice. The coverage set contains the entire mountain on 
which P lies plus parts of adjacent mountains. These points alone 
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Fig. 13 — A coverage set. 

constitute a hexagon of area 4/(7. In addition, the coverage set contains 
many smaller isolated patches on more remote mountains. These 
small patches can be so numerous that they represent most of the 
coverage area. 

If Fig. 13 represented the coverage set of a base station in a mobile 
radio telephone system, the station would only serve the hexagon of 
area 4/<r. The other small patches would lie in the service areas of 
other stations, and so these patches would represent places where the 
given station can interfere with other stations. 

As in Fig. 11, let Pi be a given peak and Qi another point at ground 
level. Suppose the distance r from Pi to Q% is known, i.e., the angle 
z = r/a in Fig. 11 is given. Let f(r) denote the probability that a 
line-of-sight path P1Q2 exists. Since an element of area dA{Qi) at Q2 
belongs to the coverage set of Pi with probability /(r)cL4.(Q 2 ), the 
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mean area covered is 

C = f ff(r)dA(Qt) = 2*0? f*f(r) sin z 2 d Zi . (36) 

Before attempting to evaluate f(r) and C, the integral (36) will be 
given a second interpretation. Now consider Q 2 at a fixed location 
and count the number of peaks visible from Q 2 . The probability that 
an element of area dA(Pi) contains a peak Pi visible from Q 2 is 
<rf(r)dA (Pi). Then the mean number of peaks visible from Q 2 is 

£ (visible peaks) = a f f f(r)dA(Pi) 

= 2ira 2 <r / / f(r) sin z 2 dz 2 

= aC. (37) 

Equation (37) can be used to derive very simple bounds on C. 
Clearly, more peaks are visible from a point Q 2 at high altitudes than 
at low. If Q 2 were itself a peak, the mean number of visible peaks 
would be #(2 ), given by (26). But Q 2 has probability zero of being 
exactly at a peak. If Q 2 is at any slightly lower altitude, Q 2 is on the 
side of a hill which obscures 180 degrees of the view from Q 2 . Thus, 
# (visible peaks) ^ iE(2 ), and (37) shows 

C g tf(2o)/(2(r). (38) 

Curiously, the right-hand side of (38) is exactly the mean area within 
the horizon as given by (35). Then (38) is a bound that was obtained 
in Section VII. 

At the other extreme, Q 2 might be on the inner sphere, where no 
peak is visible. In most cases, that event will be so unlikely that it 
will be safe to say that the worst reasonable possibility is that Q 2 is 
down in a valley near the point where three mountains meet. Here, 
the three mountain peaks are visible and so one concludes C ^ 3/cr. 

To obtain /(r), and hence C, recall that (20) is a formula for the 
path probability p(Qi, Q 2 ), depending on the altitude y at Q%. To get 
f(r) one may average p(Qi, Qi) over y (or 72). This average may be 
expressed as a sum of two terms which account for the possibilities 
that Q2 belongs to the same mountain M(P\) as Q\ or to a different 
mountain. 

In Fig. 11, if 72 = 2z, then Q 2 lies on the mountain M(Pi). This 
event has probability exp { -27nra 2 (l — cos z 2 )}. The path probability 
p(Qi, Qi) = 1 if r ^ ($t - 6)a. If r > (%w - d)a, then Q 2 lies on the 
inner sphere, the path P1Q2 is blocked, and p(Qi, Qi) = 0. 
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If 72 < 22, then Q 2 lies on a different mountain M(P 2 ). This possi- 
bility contributes a second term to f(r), 

I' 2 p(Qi,Qi)d{l - exp(-aA 2 )}, 

771=0 

where A 2 and p(Qi, Q 2 ) are given by (16) and (20). 
For r ^ (^tt — d)a, the two terms combine into 

f(r) = exp {-27nra 2 (l _ cosz 2 )l 

+ ('* (1 - ^A)exp[-a(A + A 2 )>dA 2 , (39) 

where (16) and (17) provide A 2 , A . A similar formula applies when r 
is larger, but f(r) is very small at such large ranges. 

One could find/(r) to any desired accuracy by evaluating the integral 
in (39) numerically. Instead, (39) will be replaced by a simpler approxi- 
mate formula. Since a is large, the first term of (39) is approximately 
exp ( — ffirr 2 ) ; also, A 2 = ttx 2 where x = ay%. The approximations to <p 
and A which follow are not uniformly good but are intended to apply 
in situations that contribute most of the coverage area C. Except at 
very short ranges r, a typical blocking region is more elongated than 
that shown in Fig. 11. Figure 14 is more typical. Then <p 2 = %t, 
approximately. With that approximation, the shaded region in Fig. 14 
has area A + \A 2 . It consists of a triangle, of area xr, and two extra 
lens-shaped pieces. The two extra pieces can fit together into one lens 
of exactly the shape of the blocking region for two peaks at separation 
r. Then the two extra pieces combined have area r 3 /(6ga), as in (23). 
Now the exponent <r(A + A 2 ) in (39) is approximately (r/D) 3 + \irax 2 . 

To substitute these approximations into (39), write 



Y = airr 2 

X = r 3 /(Qga) = (r/D) 3 = [r(5/3)F/^(S )] f 



(40) 




Fig. 14 — Approximation of At, and <pi. 
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and then 



/(r) = exp (— F) + xo-exp ( — X) I exp \—<rxr — %aTr 2 }rdr 

Jo 

f(r) = exp (-F) + exp (-Z){1 - exp (-FQ + i-*])} 

+ exp(-X + §F/ 1 r 2 )(2F/7r)* 

X {erf (F*[l + 7T- 1 ]) - erf (FV*-)}. 



(41) 



Figure 15 shows curves of f(r). The ordinate F* = ((nr)*r was used 
as a convenient normalized range. It may be interpreted as the square 
root of the mean number of other peaks within range r of Pi. As (40) 
shows, the parameter 2£(2 ) enters into f(r) through the variable X. 
The f(r) curves for different values of -B(S ) lie close together at short 
ranges. As the range increases, f(r) falls more sharply for small E(2o). 
There is a limiting curve, as #(2 ) — »*>, which is obtained from (41) 
by setting X = 0. As (40) shows, X = also corresponds to the limit 
a — >«> ; i.e., this limit represents the planar model. 
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Fig. 15 — Probability /(r) that a random point at ground level is visible from a peak 
r miles away. 
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Table VII — Mean coverage area C 



E(Z ) 


<jC 


25 


8.3 


50 


12.2 


100 


17.3 


200 


23.5 


400 


30.6 


800 


38.6 


1600 


47.4 



The approximations which led to (41) are poor when r is small. 
Fortunately, the planar model is good for smaller r. The curve labeled 
"exact (planar)" in Fig. 15 was obtained by a numerical integration, 
using an exact equation (39) for the planar model. The planar curve 
crosses the 0.5 probability level at Y 1 = 2.3. Then r = 1.3tf~* is the 
range at which the odds of finding a clear path are even. 

The behavior of f(r) for large r can be obtained by replacing the 
error functions in (41) by their asymptotic expansions. The leading 
terms are 

f(r) ~ exp (- Y) + **Y~ l exp (-X). (42) 

In Fig. 14, the f(r) curve starts to depart from the limiting curve at 
values of Y near J?(S ). For larger Y, the factor exp ( — X) in (42) 
becomes small rapidly. In the planar model, exp ( — X) = 1 for all Y; 
(42) then shows that f(r) ~ tt 2 /Y = r/(<rr 2 ). 

To good approximation, the integral (36) for the mean coverage 
area can be replaced by 

C = rf(r)d(m*) = a- 1 r f{r)dY. (43) 

Jo Jo 

The main contribution to C in (36) comes from the range ^ r 
= {h* ~ ") a , m which (39) holds exactly and (41) approximately. 
Then (41) will be used for f(r) in (43) and, since f(r) — * rapidly for 
larger r, the range of integration has been extended to ^ Y < » . 

Since (40) and (41) express /(r) in terms of Y and the single param- 
eter .E(2o), (43) shows that aC is a function of S(2 ) only. Table VII 
gives values of this function, obtained from (41) and (43) by numerical 
integration. These values also represent mean numbers of peaks 
visible from a random point on the ground, as (37) showed. 

In Table VII, aC appears to be a slowly increasing function of 
E(2 ). As E(2 ) — >co f the model becomes planar and then/(r) ~ ir 2 /Y. 
Then (43) shows aC — * «> in the planar model limit. 
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